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SUMMARY 


The purpose of this report is to enhance understanding of the fracture process 
of metal matrix conposites, boron/aluminum composites specifically. These composites 
experience widespread yielding of the matrix material during the fracture process. 
This plasticity coirplicates the analysis of the composites containing notches. This 
report presents experimental results for five laminate orientations of boron/aluminum 
composites containing either circular holes or crack-like slits. Specimen stress- 
strain behavior, stress at first fiber failure, and ultimate strength were deter- 
mined. Radiographs were used to monitor the fracture process. The specimens were 
analyzed with a three-dimensional elastic-plastic finite-element model. 

The first fiber failures in notched specimens with [±45]2gf [0/±45]g, and 

[02/±45]g laminate orientations occurred at or very near the specimen ultimate 
strength. For notched unidirectional specimens, the first fiber failure occurred at 
approximately one-half of the specimen ultimate strength. Acoustic emission events 
correlated with fiber breaks in unidirectional composites, but did not for other 
laminates. Circular holes and crack-like slits of the same characteristic length 
were found to produce approximately the same strength reduction. The predicted 
stress-strain responses and stress at first fiber failure compared very well with 
test data for laminates containing 0® fibers and reasonably well for [±45] 2 g 
laminates. 


INTRODUCTION 

Metal matrix composites have several inherent properties, such as high 
stiff ness-to-weight ratios and high strength-to-weight ratios, which make them 
attractive for structural applications. These composites also have greater trans- 
verse strength, a higher operating temperature range, and better environmental resis- 
tance than current epoxy-resin-matrix composites. Like polymer matrix composites, 
however, metal matrix conposites are very notch sensitive. The degree of this sensi- 
tivity depends on the notch size and shape, as well as the laminate orientation. 
Unlike typical polymer matrix composites, metal matrix composites may exhibit wide- 
spread yielding of the matrix before laminate failure. To design damage -tolerant 
structures (or to simply understand the effects of fastener holes), the laminate 
fracture strength must be known for a wide range of ply orientations, notch geom- 
etries, and loading conditions. Furthermore, a method for predicting fracture 
strength is needed to avoid testing all laminate, notch, and loading combinations of 
interest. 

The purpose of this research was to investigate the fracture process of 
boron/aluminum (B/Al) composites containing circular holes and crack-like slits. 

B/Al composites with slits had already been tested and fracture data published 
(ref. 1). To further understanding of the B/Al fracture process, four laminate 
orientations of B/Al conposites with circular holes were tested during the present 
investigation. One unidirectional monolayer specimen with a center slit was also 
tested. Stress and overall strain data were obtained from initial loading to laimi- 
nate failure, and the specimens were radiographed periodically during loading to 
record damage initiation and propagation. 



A three-dimensional elastic-plastic finite-element analysis program was modified 
to analyze the mechanical behavior of the specimens. Overall stress-strain responses 
and first fiber failures were predicted for the notched specimens. To evaluate the 
analysis, the elastic-plastic predictions were compared with experimental data. 


SYMBOLS 


AE 

acoustic emission 


a 

half-length of slit or radius of circular hole. 

S 

laminate stress, MPa 


B 

> 

> 

voliime fractions of fibers and 

matrix 

w 

specimen width, mm 



Cartesian coordinates (i = 1, 2 

, 3) 

•H 

IX 

material coordinates (i = 1, 2, 

3) 

e 

constituent strain tensor, m/m 


e 

overall laminate strain tensor. 

m/m 

a 

constituent stress tensor, MPa 


o 

overall laminate strain tensor. 

MPa 

Subscripts: 



f 

fiber 


m 

matrix 


ult 

ultimate 


1,2 

transverse directions 


3 

0° fiber, or axial, direction 



EXPERIMENTS 

A considerable amount of experimental fracture data have been generated during 
the past decade. References 2 to 4 are good examples. The past tests focused pri- 
marily on the development of stress-strain response data and laminate ultimate 
strength. The process of fracture (i.e., when damage first initiates and how it 
grows under increasing static load) has not been investigated. Therefore, the pur- 
pose of the experimental tests presented in this report is to characterize the frac- 
ture process. The overall stress-strain response and laminate failure load were, of 
course, recorded. 
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Material and Specimens 


The B/Al laminates were made by diffusion-bonding sheets of 0 . 142-mm-diameter 
boron fibers and 6061 alumin\am foil. The laminates were used in the ”as fabricated” 
condition; that is, they were not heat treated- The laminate orientations were [0]^/ 
[0/±45]g, and [±45]2g* The fiber volume fractions were 0.50 for the [0]^ 
laminate and 0-45 for the other laminates- Fracture and tensile specimens were cut 
from each sheet of material- The tensile specimens had been previously tested by 
Sova and Poe (ref- 5) to measure tensile stress-strain behavior for unnotched lami- 
nates- Fracture specimens containing crack-like slits had also been tested previ- 
ously by Poe and Sova (ref. 1). 

The fracture specimens tested in this study were 51 and 102 mm wide. Circular 
holes were machined at the center of each fracture specimen with an ultrasonic rotary 
diamond- impregnated- CO re drill. The diameter of the holes ranged from 3. 18 to 
25.4 mm. Table I lists hole diameters for each specimen width. The sketch of the 
fracture specimen in figure 1 identifies fiber orientation angles and the location 
of the linearly varying displacement transducer (LVDT) used to measure specimen 
elongation. 

One specimen consisting of a monolayer of unidirectional B/Al was also tested- 
The monolayer, 76 mm wide, contained a center crack-like slit 19 mm long and had a 
fiber volume fraction of 0.30. 


Test Procedure and Equipment 

The specimens were tested in a hydraulically actuated, closed-loop, servo- 
controlled testing machine. The load, measured by a conventional load cell, was used 
for the feedback signal. In addition, some tests were strain controlled. The loads 
(or strains) were programmed by an analog ramp generator to vary linearly with time. 
The applied load and the overall strain as measured from the LVDT were recorded on an 
X-Y plotter. The overall strain is the displacement measured by the LVDT (fig. 1) 
divided by 89 mm. 

During several tests, specimens were radiographed with an industrial-type "soft" 
X-ray machine with a 0 . 25-mm-thick beryllium window and a tungsten target. The volt- 
age was set at 50 kV and the current at 20 mA. The window of the X-ray tube was 
300 mm from the specimen. A Kodak high resolution photographic plate (type 1A) was 
mounted to the specimen on the opposite side. Exposure time was about 5 minutes. 

The load (or strain) was held constant during each radiograph. 

Specimens were radiographed either at several increments of load or when damage 
growth was indicated by acoustic emission ( AE) signals. Acoustic emission signals 
were also recorded during these tests. A Dunegan/Endevco 3000 series AE system moni- 
tored the AE events that had peak amplitudes above 65 dB. These AE events were com- 
pared with the radiographs to correlate the fiber breaks with the acoustic emissions. 


Test Results and Discussions 

Monolayer test .- Figure 2 shows a series of radiographs taken at increasing load 
levels applied to the monolayer test specimen. Only the left end of the center slit 
is shown. The dark, thin lines in the radiographs are the tungsten cores of the 
boron fibers. The light area around the tungsten core is the boron. The gray area 
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between is the aluminum matrix. Obviously, the fibers are not \miformly spaced. As 
the specimen was slowly loaded, the number of AE events in the range from 65 to above 
100 dB were recorded for each radiograph taken. ^ The number of events above 100 dB 
was compared with the number of fibers broken. 

Figure 2(a) shows the specimen at a stress level of 196 MPa with no apparent 
fiber failure nor any AE events above 100 dB. At a stress level of 216 MPa, a single 
event above 100 dB was recorded and one broken fiber was observed at the right end of 
the slit. Figure 2(b) shows two broken fibers at the notch tip at a stress level of 
237 MPa. Two AE events above 100 dB were recorded. Figure 2(c) shows six broken 
fibers with three AE events above 100 dB. Acoustic emission equipment may not record 
the sound of every fiber break, especially if breaks occur nearly simultaneously. 
Figure 2(c) clearly shows that the boron fibers break before matrix failure. This is 
because the fibers have a 0.6 percent strain to failure, whereas the matrix material 
has more than 7.5 percent strain to failure. Figure 2(d) shows the damage state just 
before failure; there were 27 broken fibers and 1 1 AE events above 100 dB. 

The results of the test of the monolayer specimen indicated a qualitative corre- 
lation between the AE events above 100 dB and the 0® fiber breaks for unidirectional 
laminates. The monolayer test illustrated that the boron fibers fail well in advance 
of the matrix material. Fiber failure first occurred at a stress of approximately 
60 percent of the ultimate load. The crack grew in a self-similar manner until the 
specimen failed. 

The [O ]0 tests.- Nine laminate specimens with circular holes were loaded 

to failure. The gross section ultimate stress for each specimen plotted versus the 
ratio of hole diameter to specimen width (2a/w) is shown in figure 3. The two speci- 
men widths are distinguished by open and closed symbols. Also plotted in figure 3 
are the ultimate strength data obtained from similar specimens containing crack-like 
slits (ref. 1), for which 2a/w is the ratio of slit length to specimen width. The 
curve for the net section ultimate strength was found by multiplying the unnotched 
tensile strength of the laminate by ( 1 - 2a/w) . This curve represents the specimen 
gross stress that corresponds to a net section stress equal to the unnotched tensile 
strength of the laminate. The unnotched tensile strengths were determined by Sova 
and Poe (ref. 5) and are listed in table II for each layup. The strength data and 
the curve for net section ultimate strength differ because of the notch effect. Note 
that at a given value of 2 a/w, the ultimate stresses are nearly the same for speci- 
mens with holes and with slits. As a result, the notch effect appears to be a func- 
tion of notch length rather than the notch radius or width. This may be attributed 
to the widespread yielding that occurs near the notch long before the failing load is 
reached. Examples of the yield patterns are presented in the analytical section. 

The wider specimens failed at a lower laminate stress for the saime value of 2a/w. 

Figure 4 shows a series of four radiographs of the [0]^ laminate and the associ- 
ated number of AE events above 100 dB. The arrow in figure 4(a) indicates the region 
of developing damage, shown as light areas indicating breaks in the fibers near the 
edge of the hole. This damage was more evident in the original radiographs. Other 
apparent fiber damage around the hole was caused by machining and was present at the 
beginning of the test. A radiograph was taken at a slightly lower load than fig- 


^The dB levels are a function of the basic AE unit and transducer used. A 
Dunegan/Endevco S 750B transducer was used with a model 1801 preamplifier which 
contained a 450B band-pass filter. 
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ure 4(a). There were no AE events and no evidence of fiber failure. The first AE 
event above 100 dB corresponded to the first observed fiber damage. The number of 
fiber breaks increased with additional AE events. Again, the AE events gave a good 
qualitative indication of fiber damage. The first fiber failure observed was at 
approximately 40 percent of the ultimate laminate stress in the location indicated by 
the arrow. As in the monolayer specimen, the damage propagated in a self -similar 
manner until failure, as can be seen in figures 4(a) through 4(d). 

The [±45] 2 s tests.- Eight [+45] laminate specimens with circular holes were 
loaded to failure. ®ie ultimate strength of each specimen is plotted versus 2a/w 
in figure 5. The strengths of specimens containing crack-like slits (ref. 1) are 
also plotted. There does not appear to be a close correlation between slit and hole 
data as was observed in figure 3 for the [0]^ laminates. This lack of correlation 
may be attributed to relative scatter in the strength data, possibly due to wide- 
spread matrix yielding. The actual range of scatter (i.e. , 100 MPa) is about the 
Scime as for the [0] ^ specimens; however, the absolute value of strength is much lower 
for the [+45] 2 g specimens. 

The unnotched tensile strength of the [±45] 2 g laminates was very width depen- 
dent, as shown in table II. No width dependency was found for the unnotched tensile 
strength of the other layups tested (ref. 1). Two net section strength lines are 
shown in figure 5 for 102-mm-wide and 51-mm-wide specimens. 

Radiographs were taken of several specimens during testing. Figure 6 shows the 
typical failing sequence during loading and the accompanying AE events. The specimen 
shown was tested under strain control. First fiber failure was not observed until 
the laminate was loaded to 95 percent of its ultimate strength, as shown in fig- 
ure 6(a). Notice that 350 AE events above 100 dB had accumulated at this stress 
level, yet only two broken fibers could be seen. AE events above 100 dB and broken 
fibers do not seem to correlate for the [±45] 2 g laminates. Perhaps these AE events 
were caused by fiber-matrix separation. Figure 6(c) shows the critical damage state 
at maximum stress level. If load control had been used, specimen fracture would have 
occurred at the damage state shown in figure 6(c). Figure 6(d) shows unstable damage 
growth occurring at a 45® angle. 

The [0/+45]g tests.- Eight [0/+45]^ laminate specimens with circular holes were 
loaded to failure. The ultimate laminate stress for each specimen is plotted versus 
2a/w in figure 7. The hole data again agree very closely with the slit data of Poe 
and Sova (ref. 1). Like the [0] ^ laminates, the wider specimens failed at a lower 
laminate stress for the same value of 2a/w. 

Radiographs of specimens were taken at increasing load levels approaching the 
ultimate strength. The radiographs showed no signs of fiber failure prior to lami- 
nate failure. Presumably, the fibers all failed in rapid sequence at the ultimate 
load. Acoustic emission events were recorded at load levels above 25 percent of the 
ultimate load, yet no fiber failures were observed. 

The [02/±45]g tests .- Nine [02/±45]g laminate specimens with circular holes were 
loaded to failur*^ Figure 8 presents the ultimate laminate strength plotted versus 
2a/w for each specimen. As previously noted for the [0]^ and [0/±45]g laminates, 
the [02/±45]g laminates with holes and with slits exhibited virtually no difference 
in ultimate strengths, but notch length did have an effect. 

Radiographs were also taken of a [02/±45]g laminate near the ultimate stress. 

As for the [0/±45]g laminates, no fiber failures were observed before leunciinate fail- 
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ure. Again, many AE events above 100 dB were recorded throughout the tests, but none 
seemed to correlate with fiber failures. 


Summary of Experimental Data 

In summary, the experimental phase of the investigation indicated the 
following: 

1. First fiber failure occurred very near the ultimate notched Icuninate strength 

(within 5 percent) for [±45]2g/ [0/±45]g, and [02/±45]g B/Al laminates containing 

holes. Unidirectional laminates with holes e3q>erienced first fiber failures at 
approximately 50 percent of ultimate laminate stress. Poe and Sova (ref. 1) and 
Jones (ref. 6) found similar results for B/Al laminates containing crack-like slits. 

2. Acoustic emission events correlated with fiber breaks in unidirectional lami- 
nates, but did not for other laminate orientations. 

3. The ultimate strength decreased with increased specimen notch length (slit 
length or hole diameter). 


ANALYTICAL MODELING OF FIBROUS COMPOSITES 

In aluminiam matrix composites the matrix yields, whereas the fibers generally 
behave elastically until they break. Therefore, to model the behavior of metal 
matrix composites, the elastic-plastic behavior of the matrix and the elastic behav- 
ior of the fiber must be accounted for. An analysis that models this two-phase 
behavior with any specimen configuration and incorporates fiber and matrix failure 
was developed in the present study. The analysis used is a three-dimensional finite- 
element program (ref. 7) called PAFAC (plastic and failure analysis of composites), 
which was developed from a program written by Bahei-El-Din et al. (refs. 8 and 9) . 
PAFAC uses an eight-noded hexahedral element and the least -squares method to compute 
the element strain coitponents from the total of 28 strain-displacement relationships 
associated with all possible combinations of the eight vertices taken two at a 
time. Each hexahedral element represents a unidirectional conposite material whose 
fibers are arbitrarily oriented in the structural (Cartesian) coordinate system. 


Finite-Element Model 

Two typical finite-element meshes used in the PAFAC program are shown in fig- 
ure 9. The mesh used for a specimen with a center slit (fig. 9(a)) consisted of 
about 350 eight-noded hexahedral elements with about 1500 degrees of freedom. The 
mesh used to model a specimen with a circular hole in the center (fig. 9(b) ) con- 
sisted of about 400 eight-noded hexahedral elements with about 2000 degrees of 
freedom. In all cases, only one-eighth of the specimen was modeled because of 
symmetry. Symmetry conditions were not satisfied when only one-eighth of the 
laminates with ±45° plies were modeled; however, we felt that the error introduced 
by assuming symmetry for these cases was insignificant. 

At least one layer of elements is needed for every ply angle in the laminate. 
For example, a [0]^ laminate could be modeled by one layer of elements, while a 
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[±45] 2 g laminate would require at least two layers of elements. A uniform stress was 
applied to the end of the specimen to simulate the loading- 

To evaluate the accuracy of the model used, the finite-element results obtained 
using the mesh shown in figure 9{b) for an isotropic material were compared with the 
exact solution for a hole in an isotropic plate corrected for finite width. For the 
axial stress at the edge of the hole, the two solutions differed by less than 
4 percent. 


Material Model 

The PAFAC program uses a continuum material model developed by Bahei-El-Din and 
Dvorak (refs. 9 through 11) to represent the essential aspects of the elastic-plastic 
behavior of composite laminates. The material model is briefly described in this 
section. 

The model consists of an elastic-plastic matrix unidirectionally reinforced by 
continuous elastic fibers. Both constituents are assumed to be homogeneous and iso- 
tropic. The fibers are assumed to have a very small diameter, so that although the 
fibers occupy a finite volume fraction of the composite, they do not interfere with 
matrix deformation in the two transverse directions, but only in the axial (fiber) 
direction. Figure 10(a) shows a schematic of this lamina model. It can also be 
represented by parallel fiber and matrix bars or plates with axial coupling, as 
illustrated in figure 10(b). In figure 10, the fiber (axial) direction is parallel 
to the X 2 ~a 5 cis, and the x^- and x^-axes represent the transverse directions. 

If the Cartesian coordinates are chosen so that X 2 coincides with the fiber 
direction, the second-order tensors of the independent stress and strain components, 

0 and e, are expressed as 


= [c 


11 22 33 12 "13 23 


f 


e 



^22 ^33 ^12 ^13 ^23^ 


where Yij = (i,j=1, 2, 3; i^j) are the engineering shear strain 

components . 


For equilibrium and conpatibility , several requirements are imposed on the mate- 
rial model shown in figure 10. The stress average in each constituent can be related 
to the overall composite stress a in the axial (fiber) direction as follows: 
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A bar over a symbol indicates overall composite stress or strain, and the sub- 
scripts f and m denote quantities related to the fiber and matrix. The volume 
fractions v^ and v^ are such that “*■ other stress components 

in each constituent were assumed to be uniform and to obey the following equilibrium 
equations: 


^11 = 

^22 " ^'^22^f " ^‘^22^m 

®12 "" ^ ^®12V \ 

^13 ^ 

^23 "" ^*^23^f “ ^*^23^tn 


( 2 ) 


The only constraint in the model is in the axial (fiber) direction; the matrix 
and fiber must deform equally. Thus, 


^33 “ ^^33^f “ ^^33^iti 


(3) 


The other strain coitponents can be related to the overall strain e as follows: 


^ij 




^m^^ij ^m 


(ij ^ 33) 


(4) 


Since the fibers are elastic up to failure, the inelastic strains of the lamina 
are caused by matrix deformation. Because the fiber imposes an elastic constraint on 
the matrix which affects the shape of the lamina yield surface, additional kinematic 
components appear in the hardening rule of the lamina and influence the magnitude of 
the overall plastic strains. All aspects of the yield behavior were examined and 
accounted for in the formulation of the lamina constitutive equations. These equa- 
tions are explicitly described in references 9 and 11. 


Elastic-Plastic Response 

The material model described above was used to characterize the elastic-plastic 
behavior of simple composite laminates. The stress-strain curve of the matrix 
material was modeled with a Ramberg-Osgood equation (ref. 12). In addition, the 
program PAFAC is able to calculate the fiber stress in each element and to model 
fiber breakage. 

The yield pattern of the matrix material can be determined with the PAFAC 
model. Figure 11 shows the extent of yielded elements for increasing levels of 
stress applied to a [0]^ laminate with a 25. 4-mm- diameter hole and 2a/w = 0.25. 

Only one-quarter of the specimen is shown in the figure. The sketches show the 
growth of yielding as applied stress is increased. The first yielding occurred at 
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a stress of 34 MPa, shown in figure 11(a)* At the applied stress of 125 MPa, most of 
the specimen had yielded; only a small portion of the model above the hole remained 
elastic. Although not shown in the figure, the entire specimen had yielded at a 
stress of 275 MPa. For this example, the test results showed that fibers began to 
fail at about 260 MPa and the specimen failed at a stress of 540 MPa. Widespread 
yielding occurred before fibers broke or the specimen failed. This yielding influ- 
enced the ultimate fracture behavior of the composite. This example emphasizes the 
necessity for a model that can predict the elastic-plastic nonlinear behavior of the 
matrix phase and the elastic linear behavior of the fiber phase for metal matrix 
composites. 

This technique can be applied to each layer of laminates with plies at different 
angles. Figure 12 shows the yield pattern for a [0/+45]g laminate with a 25.4-mm- 
diameter hole and 2a/w = 0.25. The sketches show the extent of the yielding in the 
0® and 45® layers as the applied stress is increased. The yielding in the -45® layer 
is very similar to the pattern shown for the 45® layer. The +45® layers yielded much 
earlier in the loading history than the 0® layer. Although not shown in the figure, 
the entire specimen (all layers) had yielded at 150 MPa. As in the previous example, 
widespread yielding occurred before fibers broke or the specimen failed. 


ANALYSIS OF FAILURE 

The following equation was incorporated into the PAFAC program as the fiber 
failure criterion: 



In this equation, represents the current tensile stress in the fiber and 

represents the current shear stress in the fiber. The ultimate strengths of the 
fiber in the tensile and shear modes are <^ 33 ^ ^^3^' respectively. For the 

meshes used, the element at the notch tip was approximately three fibers wide. A 
mesh fine enough to model a fiber width by one element width would be prohibitively 
expensive. Therefore, an averaging technique was used to determine the ultimate 
strengths, ^ 33 ^ ^^^ 3 ^ fibers within one element. 

Since these two stengths were unknown, two test cases were required. The 
monolayer with a slit and a [+45] 2 g specimen containing a hole were selected. These 
two layups were chosen because for the monolayer the dominant fiber stress at the 
crack tip was the tensile stress ^^ 3 ^ while in the [+45] 2 g laminate the dominant 

fiber stress was the shear stress cr^ 3 * From the radiographs, the laminate stress at 
first fiber failure could be closely estimated. The finite-element model for each 
laminate configuration was loaded to the approximate stress when first fiber failure 
occurred. The corresponding fiber stresses 0^3 and were recorded for the 

element where first fiber failure was observed. The pairs of values for 
0^2 both the [±45] 2 g and monolayer specimens were substituted in equation (5). 

The two equations were then solved for the two unknowns, ^^3 * This gave 

the values of 2300 MPa and 270 MPa, respectively. First fiber failure was predicted 
for the other laminates using these values. 
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To validate this method of determining element fiber strengths, these values 
were used in the f inite-element model to study the fiber stress in the region of 
notches. For the unidirectional monolayer with a center slit, figure 13 shows how 
the fiber stress varies in the region of the crack tip. The solid line represents 
data results from the finite-element analysis. The program predicts stresses only 
at element centroids; therefore, to obtain the stress in the first fiber, it was 
necessary to extrapolate the stress from the element centroid to the edge of the 
notch. This extrapolation is represented by the dashed line in figure 13. As indi- 
cated in the figure, the crack-tip element contained approximately 2.3 fibers. The 
ultimate strength of a single fiber was assumed to be 3160 MPa, which was obtained 
from tensile tests of unnotched 0® B/Al specimens (ref. 5). The stress at the center 
of the fiber next to the crack tip was 3160 MPa when the applied remote stress was 
262.5 MPa, a reasonable estimate of the remote stress at first fiber failure. The 
corresponding stress at the center of the element next to the crack tip was found to 
be about 2300 MPa. Figure 14 shows a similar plot of the fiber stress concentration 
for a [0]g laminate with a hole. In this example, the element next to the hole con- 
tained approximately 3.6 fibers. Here the stress at the center of the fiber next to 
the hole was found to be about 3160 MPa when the applied remote stress was 465 MPa. 
Again, the stress at the center of the element next to the hole was about 2300 MPa. 

At each incremental load level, the program checked the state of fiber stresses 
in each element. When equation (5) was satisfied, the fibers of that element were 
broken. In the PAFAC program, it is possible to continue the analysis after first 
fiber failure. However, because of the inability of the program to model element 
unloading when a fiber broke, we felt that for this study, any calculations made 
after first fiber failure would be unrealistic; therefore, the analysis was termi- 
nated at the point of first fiber failure. 


COMPARISON OF PREDICTED AND EXPERIMENTAL RESULTS 

In this section, the predicted mechanical behavior of notched B/Al composite 
specimens is compared with actual test data. To show the capability of the analysis 
to predict small-scale, localized mechanical behavior, crack-opening displacements 
(COD) were predicted for the specimens containing slits tested by Awerbuch and Hahn 
(ref. 4). The data of reference 4, which were the only COD data available for com- 
parison, were for unidirectional specimens. A plot of the specimen applied load 
versus the measured and predicted COD is shown in figure 15. The PAFAC model pre- 
dicted COD very accurately up to first fiber failure. 

As discussed earlier, the finite-element program PAFAC can predict first fiber 
failure around a notch. The first fiber is assumed to fail when the fiber stresses 
within an element exceed the criterion given in equation (5). As previously 
explained, the elements near the notch tip were approximately three fibers wide. 
Therefore, "first fiber failure" actually indicated failure of all three fibers 
within an element. The following sections show the correlation between the pre- 
dicted first fiber failure and the observed first fiber failure for several dif- 
ferent laminates containing both holes and slits. Ultimate strength data are also 
presented. 


Specimens With Holes 

In figure 16 the applied laminate stress is plotted against the specimen overall 
strain for a unidirectional specimen containing a circular hole. The solid line 
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represents the predicted response up to first fiber failure. There is excellent 
correlation between the predicted and observed results up to first fiber failure- As 
previously noted, the first fiber failed quite early in the loading history for the 
unidirectional specimens. The sketch shows the predicted location of the first fiber 
failure. 

Figure 17 shows the predicted curve and ejqjerimental results for two [±45] 2 ^ 
specimens. There was considerable scatter in the data. The predicted curve lies 
just above the upper e3q)erimental curve, for which the stress at first fiber failure 
is expected to be near the laminate failure stress, as observed in the lower curve. 

The results for a [0/±45] ^ laminate are shown in figure 18. Again, there is a 
very good correlation between the predicted and the experimental stress-strain 
behavior. First fiber failure was predicted to occur near laminate failure- Recall 
that first fiber failure was not observed for this laminate and therefore was assumed 
to occur at laminate failure. The first fiber failure was predicted next to the hole 
in the 0® ply element on the transverse axis, as shown by the sketch. Figure 19 
shows similar results for a [02/±45] ^ laminate- This layup was more dominated by 
0® fibers than the [0/±45] ^ laminate, so less plasticity was predicted. 


Specimens With Slits 

The PAFAC analysis was also used to predict first fiber failure of a unidi- 
rectional monolayer containing a crack-like slit. The results are shown in fig- 
ure 20. The yield pattern that developed at the notch tip (fig. 20(a)) extends 
perpendicular to the notch, parallel to the fibers. This behavior has been observed 
by Jones (ref. 6). Figure 20(b) presents a plot of the normalized fiber stress 
(fiber failure criterion) versus laminate stress. The solid line represents the 
predicted fiber stress up to first fiber failure. The dashed line shows the laminate 
stress for which the first fiber failure was observed. Again, the predicted and 
observed laminate stresses at first fiber failure agree very well. 

Several specimens containing crack-like slits tested by Poe and Sova (ref. 1) 
were also analyzed with the PAFAC program- Figure 21 shows the predicted and experi- 
mental results for a 51-mm-wide [0]^ specimen containing a 15-mm slit. Poe and Sova 
measured the remote strain with strain gages- The predicted strains at the same 
location are compared in figure 21 with the measured remote strains- Again, as with 
all the other unidirectional specimens tested, the first fiber is predicted to fail 
early in the loading history. Poe and Sova observed this early failure of fibers but 
did not record the stress at which the fibers began to fail. The stress-strain pre- 
diction is quite good. 

Figure 22 shows similar analytical and experimental results (ref. 1) for a 
102-mm-wide [0/±45]g specimen containing a 10-mm slit. The PAFAC model predicted 
stress-strain response very well. The predicted first fiber failure was very close 
to the laminate failure stress. 


Summary of Analytical Predictions 

The PAFAC model produces very good predictions of observed mechanical behavior 
of B/Al conposite specimens containing either circular holes or crack-like slits. 
The model accurately predicts stress-strain responses of the notched composites, 
including widespread yielding, and can predict first fiber failure. Because first 
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fiber failure occurs very close to overall laminate failure in all laminates examined 
except the unidirectional ones, the first fiber failure criterion is sufficient for 
predicting ultimate strength of most laminates containing notches. Figures 23 and 24 
summarize the elastic-plastic strength predictions made with the first fiber failure 
criterion. The predictions were made on specimens of various widths and ratios of 
notch length to specimen width (2a/w), as noted in the figures. 


CONCLUSIONS 

Metal matrix composites in general, boron/aluminum composites in particular, are 
subject to widespread yielding. T^is plasticity conplicates the analysis of metal 
matrix conposites containing notches. This report presented experimental results for 
five laminate orientations of boron/aluminum conposites containing circular holes and 
crack-like slits. %>ecimen stress-strain behavior, stress at first fiber failure, 
and ultimate strength were determined from tests. The specimens were then analyzed 
with a three-dimensional elastic-plastic finite-element program. Predicted mechani- 
cal behavior was conpared with test results. The following conclusions were derived 
from this work: 

1. For notched [±45]2g/ [0/±45]g, and [ 02 /^^^^s specimens, the first fibers 

failed at a stress level near the specimen ultimate strength. 

2. For notched unidirectional laminate specimens, the first fiber failure 

occurred at approximately one-half of the specimen ultimate strength. 

3 . Acoustic emissions were useful for monitoring fiber breaks in unidirectional 

composites, but were inconclusive for other laminate orientations. 

4. Circular holes and crack-like slits of the same characteristic length pro- 

duced approximately the same strength reduction. 

5. Matrix plasticity had to be accounted for to accurately predict the stress- 

strain and failure responses of boron/aluminum composites. 

6. The PAFAC model predicted stress-strain responses and stress at first fiber 

failure very well for laminates containing 0® plies and reasonably well for 

[+45] 2 g laminates. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
July 8, 1983 
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Figure 1.- Typical B/Al composite tensile specimen with linearly varying 
displacement transducer ( LVDT) attached. 







(a) No damage - 


Figure 2.- Fracture process of a unidirectional monolayer containing a crack- 
like slit. = 0.30; w = 78 mm; 2a = 19 mm. (Figure shows one end of 

the center slit. ) 
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(b) Two broken fibers. 
Figure 2.- Continued. 
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(c) six broken fibers and matrix crack. 
Figure 2.- Continued. 
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(d) Twenty-seven broken fibers and matrix crack. 
Figure 2.- Concluded. 
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Figure 3.- Strength of [0]g laminates with holes or slits 
V- = 0.50. (Slit data from ref. 1.) 
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(a) First evidence of damage at 270 MPa. 

Figure 4.- Fracture process of a [0]^ laminate containing a circular hole. 

= 0.50; w = 102 mm; 2a = 12.7 mm. 
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(b) Damage at 400 MPa- 


Figure 4-- Continued. 
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(c) Damage at 550 MPa. 
Figure 4-- Continued - 
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Figure 5 .- Strength of [± 45 ] 2g laminates with holes or slits. 

( Slit data from ref . 1 . ) 









' 1.0 mm' 


(a) First evidence of damage at 120 MPa. 


Figure 6.- Fracture process of a [±45] laminate containing a circular 
hole. = 0.45; w = 51 mm; 2a = 12.7 mm. (Test was conducted 

under strain control. ) 
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(b) Damage at 122 MPa. 
Figure 6.- Continued. 
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(c) Unstable damage state- 
Figure 6.- Continued- 
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(d) Unstable damage growth. 
Figure 6.- Concluded. 
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Figure 7.- Strength of [0/±45]g laminates with holes or 
slits- = 0.45. (Slit data- from ref. 1.) 
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(a) Slit. 


(b) Hole- 


Figure 9.- Typical meshes used in finite-element analysis program. 



(a) Schematic model of lamina- 


(b) Parallel fiber and matrix bars 


Figure 10.- Element material model. 
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(c) S = 53.33 MPa. 



(d) S = 66.67 MPa. 


Figure 


12.- Yield pattern for 
2a = 25.4 mm; w 


[0/±45]g laminate with a circular hole. 
= 102 mm; ”^61 MPa. 
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Figure 14.- Predicted fiber stress distribution near 
circular hole for [O]^ laminate* 2a = 12.7 mm; 
w = 102 mm. 
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Figure 15.- Comparison of measured and predicted crack- 
opening displacements for 1 .32-mm-thick [0]g laminate. 
2a = 7.73 mm; w = 25.4 mm. (Data from ref. 4.) 
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Figure 16.- Stress-strain response to failure of [0]g laminate 
with circular hole. 2a = 12.7 mm; w = 102 mm. 
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Figure 17.- Stress-strain response to failure of two [±45] 2g laminates 
with circular holes. 2a = 25.4 mm; w = 102 mm. 
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Figure 18.- Stress-strain response to failure of [0/±45]g laminate 
with circular hole. 2a = 25.4 mm; w = 102 mm. 
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Figure 19.- Stress-strain response to failure of [02/±45]g laminate 
with circular hole. 2a = 25.4 nun; w = 102 mm. 
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(b) First fiber failure* 


Figure 20*- First fiber failure of monolayer specimen. 
2a = 19 Ttun; w = 78 mm. 
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REMOTE STRAIN, m/m 

Figure 2U- Stress-strain response to failure of [0]^ laminate with 
crack-like slit. 2a = 15 mm; w = 51 mm. 
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REMOTE STRAIN, m/m 

Figure 22.- Stress-strain response to failure of [0/±45]g laminate 
with crack-like slit. 2a = 10 mm; w = 102 mm. 
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Figure 23.- First fiber failure compared with laminate ultimate 
stress for specimens with circular holes. w = 102 mm. 
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Figure 24.- First fiber failure compared with laminate ultimate 
stress for specimens with crack-like slits. 
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